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ABSTRACT 

We study the model-building conditions under which a sizeable Oz//3/3-decay signal to the recently 
reported level of 0.4 eV is due to Kaluza-Klein singlet neutrinos in theories with large extra dimensions. 
Our analysis is based on 5-dimensional singlet-neutrino models compactified on an S 1 /Z2 orbifold, 
where the Standard-Model fields are localized on a 3-brane. We show that a successful interpretation of 
a positive signal within the above minimal 5-dimensional framework would require a non-vanishing shift 
of the 3-brane from the orbifold fixed points by an amount smaller than the typical scale (100 MeV) -1 
characterizing the Fermi nuclear momentum. The resulting 5-dimensional models predict a sizeable 
effective Majorana-neutrino mass that could be several orders of magnitude larger than the light 
neutrino masses. Most interestingly, the brane-shifted models with only one bulk sterile neutrino also 
predict novel trigonometric textures leading to mass scenarios with hierarchical active neutrinos and 
large v^-v T and v e -Vn mixings that can fully explain the current atmospheric and solar neutrino data. 
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1 Introduction 



Recently, realizations of phenomenologically viable theories with large compact dimensions 
of TeV size j]Q have enriched dramatically our perspectives in searching for physics beyond the 
Standard Model (SM). Among the possible higher-dimensional realizations, sterile neutrinos 
propagating in large extra dimensions OEHHEllSllZllHlini may provide interesting alternatives 
for generating the observed light neutrino masses. On the other hand, detailed experimental 
studies of neutrino properties may even shed light on the geometry and/or shape of the new 
dimensions. In this context, one of the most sensitive experimental approaches to neutrino 
masses and their properties is the search for neutrinoless double beta decay [10]. Neutrinoless 
double beta decay, denoted in short as QuflP, corresponds to two single beta decays [TT]IT2] 
occurring simultaneously in one nucleus, thereby converting a nucleus (Z, A) into a nucleus 
(Z + 2,A), i.e. 

This process violates lepton number by two units and hence its observation would signal physics 
beyond the SM. To a very good approximation, the half life for a Ou/3/3 decay mediated by light 
neutrinos is given by 

K?r i = ^iawg 01 , (i.i) 

where (m) denotes the effective neutrino Majorana mass, m e is the electron mass and -Mov/3/3 
and Goi denote the appropriate nuclear matrix element and the phase space factor, respectively. 
For details, see [TU1III1II2] and our discussion in Section 4. 

Most recently, the Heidelberg-Moscow collaboration has reanalyzed its experimental 
data J3]> using appropriate statistical methods as well as new information from the form 
of the contributing background. They found an excess of Ov/3/3 events, with statistical signifi- 
cance 2.2-3.1 a depending on the method used. From this result, a half-life of 1.5io 6 y 8 x 10 25 
years at 95% confidence level (CL) for 76 Ge is deduced, which implies an absolute value for the 
effective Majorana- neutrino mass: 

|<m) | = 0.39lg;t|eV (95% CL) , (1.2) 

allowing an uncertainty of the nuclear matrix element values of ±50%. 

The above experimental result (jl.2|) . combined with information from solar and atmo- 
spheric neutrino data, restricts the admissible forms of the light-neutrino mass hierarchies 
in 4-dimensional models with 3 left-handed (active) neutrinos. The allowed scenarios contain 
either degenerate neutrinos or neutrinos that have an inverse mass hierarchy [Tl| . Evidently, 
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a successful interpretation of a positive 0uf3f3 signal of the appropriate size mentioned above 
imposes certain constraints on the structure of a theory. Here, we study these constraints on 
the model building of minimal 5-dimensional theories compactified on a. S 1 /Z 2 orbifold. Within 
the framework of theories with large extra dimensions, previous studies on neutrinoless double 
beta decays were performed within the context of higher-dimensional models that utilize the 
shining mechanism from a distant brane jT^j and of theories with wrapped geometric space [H)] . 
In Ref. jT3], the Qvf3(3 decay is accompanied with emission of Majorons, whereas the prediction 
in falls short by two orders of magnitude to account for the observable excess in ()1.2j) . 

In this paper, we consider an even more minimal higher- dimensional framework of lepton- 
number violation, in which only one 5-dimensional (bulk) sterile neutrino is added to the field 
content of the SM. In this minimal model, the SM fields are localized on a 4-dimensional 
Minkowski subspace, also termed 3-brane. The violation of the lepton number may oc- 
cur in three distinct ways: (i) by adding lepton-number violating bilinears of the Majorana 
type in the Lagrangian; (ii) by generating lepton-number-violating mass terms through the 
Scherk-Schwartz mechanism ; (iii) by simultaneously coupling the Z2-even and Z 2 -odd two- 
component spinors of the 5-dimensional sterile neutrino to the same left-handed charged lepton 
state. As we will see in Section 2, the last case (iii) is only possible if the 3-brane describing our 
observable world is shifted from the S 1 /Z 2 orbifold fixed point. Here, we should also note that 
after integration of the extra dimension, the 5-dimensional orbifold model predicts an infinite 
tower of Kaluza-Klein (KK) neutrinos, for which the cases (i) and (ii) become fully equivalent. 

One salient feature of the S 1 /Z 2 orbifold compactification is that the KK neutrinos group 
themselves into approximately degenerate pairs of opposite CP parities. As a result, the lepton- 
number-violating KK- neutrino effects cancel each other and so the predicted 0^/3/3 decay turns 
out to be exceedingly small to account for the recent observable excess. The latter appears to 
be a major obstacle in theories with large extra dimensions and imposes by itself constraints 
on the model-building of higher- dimensional theories. A minimal way that avoids the above 
disastrous CP-parity cancellation effects on the OvfiP decay amplitude would be to arrange 
the opposite CP-parity KK neutrinos to couple to the W ± bosons with unequal strength. 
Within the minimal 5-dimensional orbifold model outlined above, such a realization can be 
accomplished only if the 3-brane is displaced from one of the S 1 /Z 2 orbifold fixed points. In 
our phenomenological bottom-up approach, the amount of brane- shifting is not arbitrary but 
dictated by the requirement that the model can accommodate the result (jl.2j) for the effective 
Majorana-neutrino mass. In particular, we will see in Section 4 how the resulting brane-shifted 
5-dimensional models can predict a sizeable effective Majorana-neutrino mass that could be 
several orders of magnitude larger than the light neutrino masses and hence than the difference 
of their squares as required from neutrino oscillation data. 
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Another important constraint on the structure of higher- dimensional neutrino theories arises 
from their ability to explain the solar and atmospheric neutrino data by means of neutrino 
oscillations. In particular, orbifold models with one bulk neutrino, as those considered earlier 
in the literature [21IH0IBlin] , seem to prefer the Small Mixing Angle (SMA) Mikheev-Smirnov- 
Wolfenstein (MSW) solution which is highly disfavoured by recent neutrino data analyses. 
Alternatively, if all neutrino data are to be explained by oscillations of active neutrinos with 
a small admixture of sterile KK component, then the compactification scale has to be much 
higher than the brane-Dirac mass terms. After integrating out the bulk neutrino of the model, 
the effective light-neutrino mass matrix has a rather restricted form; it is effectively of rank 1. 
As a result, two out of the three active neutrinos are massless. This is rather undesirable, since 
only one neutrino-mass difference can be formed in this case, so accommodating all neutrino 
oscillation data proves rather problematic jZUHUHl- However, the earlier studies have not included 
the possibility of a shifted brane. As was mentioned above, brane-shifting gives rise to sizeable 
lepton-number violation. Hence, the tree-level rank-1 form of the effective neutrino mass matrix 
can be significantly modified through lepton-number violating Yukawa terms. As we will see in 
Section 5, the resulting neutrino mass matrix has sufficiently rich structure to enable adequate 
description of the neutrino data. 

Our paper is organized as follows: Section 2 describes the low-energy structure of the 5- 
dimensional orbifold models. Technical details are relegated to the appendices. In Section 3, 
we study the renormalization-group (RG) effects of the neutrino Yukawa couplings and their 
possible impact on the QvfiP decay amplitude. In Section 4 we give estimates of the effective 
Majorana-neutrino mass, which are predicted in these models presented in Section 2. In Section 
5, we discuss the compatibility of such models with solar and atmospheric neutrino data. 
Finally, we draw our conclusions in Section 6. 

2 Minimal higher-dimensional neutrino models 

In this section, we will describe the basic low-energy structure of minimal higher- dimensional 
extensions of the SM that include singlet neutrinos. In particular, we assume that singlet 
neutrinos being neutral under the SU(2) i ®U(l)y gauge group can freely propagate in a higher- 
dimensional space of [1 + (3 + 5)] dimensions, the so-called bulk, whereas all SM particles are 
localized in a (1 + 3)-dimensional subspace, known as 3-brane or simply brane. However, even 
singlet neutrinos themselves may live in a subspace of an even higher-dimensional space of 
[1 + (3 + n g )} dimensions, with 5 < n g , in which gravity propagates. 
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We shall restrict our study to 5-dimensional models, i.e. the case 5 = 1, where the singlet 
neutrinos are compactified on a S 1 /Z 2 orbifold. Specifically, the leptonic sector of our 5- 
dimensional model consists of the SM lepton fields: 

L(x) = W, (2-1) 

with I = e, n, t, and one 5-dimensional (bulk) singlet neutrino: 

\ vfav) J 

where y denotes the additional compact dimension, and £ and rj are 5-dimensional two- 
component spinors. The SM leptons are localized at the one of the two fixed points of the 
S 1 /Z 2 orbifold, e.g. y = 0. For generality, we will assume that the brane is shifted from the 
orbifold fixed point to y = a. 

As usual, we impose the periodic boundary condition N(x, y) = N(x, y + 2nR) with respect 
to y dimension on the singlet neutrino field. In addition, the action of S l /Z 2 orbifolding on the 
5-dimensional spinors £ and rj entails the additional identifications: 

€(x,y) = £(x,-y), r){x,y) = -r){x,-y). (2.3) 

In other words, the spinors £ and r\ are symmetric and antisymmetric under a y reflection, 
respectively. 

With the above definitions, the most generic effective 4-dimensional Lagrangian of such a 
model is given by pIJEj 1 

C cS = J dy jiV(V<9 M + l5 d y )N - l - (MN T C^- l N + h.a 



+ 5{y 



h l - h l 



L$*£ + i— L$*ri + h.c. 



+ 5(y-a)C SM \, (2.4) 



{M F ) S / 2 s {M F ) S / 2 

where $ = ia 2 ^* is the hypercharge-conjugate of the SM Higgs doublet $, with hypercharge 
F($) = 1, and £sm denotes the SM Lagrangian which is restricted on a brane at y = a 0- In 
addition, Mp is the fundamental n s -dimensional Planck scale and 5 = 1 for sterile neutrinos 
propagating in 5 dimensions. Notice that the mass term m^NN is not allowed in (|2.4j) . as 
a result of the Z 2 discrete symmetry. Finally, in writing (|2.4j) . we have used the following 
conventions: 

(-? r.i. w> 




1 Further non-covariant extensions to this model have been considered in [Sj. 
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with <t m = (I2, cr) and <t m = (I2, —tr), where 01,2,3 are the usual Pauli matrices. 



We now proceed with the compactification of the y dimension of the S 1 /Z2 orbifold model. 
Because of their symmetric and antisymmetric properties (|2.Hjl under y reflection, the two- 
component spinors £ and rj can be expanded in a Fourier series of cosine and sine harmonics: 



rj(x,y) 



1 



1 



£o(z) + 
00 

Vn(x) sin 



1 ^ , f s ny 

; > C n (X) COS 



n=l 

ny 
R 



where the chiral spinors £ n (x) and r/ n (:r) form an infinite tower of KK modes. 



(2.6) 
(2.7) 



After substituting ()2.6|) into ()2.4j) and integrating out the y coordinate, we obtain the 
effective 4-dimensional Lagrangian 



1 00 

Ces = C SM + £ o (i0^)£o + (h[ {0) L$*£ - -M£ Q £ + h.c.) + £ 



n=l 



n 



TfoTfo + h.c. 



+ V2(h l l n) L$Xn + h 2 in) L$*Vn + h.c. 



where 



l(n) 



h 1 ^ 



{2ttM f R) 5 / 2 



hi, 



(2%M F R) S / 2 



cos 



sin 



na 
~R 
na 
~R 



M F \ 5 / n 



Mv 



na 



M F \ 5 / n 



M v 



h x cos ( — 

, . , na 
ho sin — 
R 



(2.8) 

(2.9) 
(2.10) 



In deriving the last step on the RHS's of (j2.9j) and 1)2.10)1 . we have employed the basic rela- 
tion among the Planck mass Mv, the corresponding n g - dimensional Planck mass M F and the 
compactification radii R (all taken to be of equal size): 



M P = (2vr M F R) n ° /2 M F . 



(2.11) 



From ()2.9j) and (j2.10j) . we see that the reduced 4-dimensional Yukawa couplings hil can be 
suppressed by many orders of magnitude [312] if there is a large hierarchy between Mp and 
the quantum gravity scale Mp. Thus, if gravity and bulk neutrinos feel the same number of 
extra dimensions, i.e. 5 = n g , the 4-dimensional Yukawa couplings hi and J12 are naturally 
suppressed by a huge factor M F /Mv ~ 10~ 15 , for M F « 10 TeV. From 1)2.2)1 . we observe that 
£ and fj belong to the same multiplet and hence have the same lepton number. It then follows 
from ()2.8)) that the simultaneous presence of h^ and J12 in an amplitude gives rise to lepton 
number violation by two units. 
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We should note that the above large suppression factor can be also obtained in a 5- 
dimensional neutrino model (5 = 1), where gravity propagates in a 6-dimensional space with 
compactification radii Ri and R 2 of unequal size (n g = 2). In this case, one has to use the 
general toroidal compactification condition: 



(2irM F ) n °' 2 {R l R 2 ...R nq fl 2 M j 



F ■ 



[2.12) 



Note that (j2.12j) reduces to (j2.11j) if all compactification radii are equal. With the help of 
(127121) . we find for n g = 2 



a 

'<a,2 



(2tiM f R 2 ) 1/2 -^fh[ 



2 ■ 



(2.13) 



We easily see that if R 2 ~ 1/Mp, the original Yukawa couplings undergo the same large 
degree of suppression by a factor M F /Mp. 

If the brane were located at the one of the two orbifold fixed points, e.g. at y = 0, the 
operator L$*r/ would be absent as a consequence of the Z 2 discrete symmetry. However, if the 
brane is shifted by an amount o ^ 0, the above operator is no longer absent. In fact, as we 
will see in Section 4, the coexistence of the two operators and L$*r/ breaks the lepton 

number leading to observable effects in neutrinoless double beta decay experiments. 



Let us now introduce the weak basis for the KK-Weyl spinors 

X±n = -j= ( in ± Vn ) , 

which enables to express the effective kinetic term of the neutrino sector as follows: 

1 



kin 



where x 1 



(yi, fo, xij x-ij 
/ o 



KK 



m 
m M 



rn 




rn 

m 

m 

m 

V : : 



M+i 









2 

. . ) and 

m 



M-i 





(2.14) 



(2.15) 



m 






m 





M — 



(2.16) 



J 



with m = vhi/y/2. In a three-generation model, m and /ii are both 3-vectors in the flavour 
space, i.e. h% = (h\, h^, h\) T ■ We will discuss intergenerational mixing effects in more detail 
in Section 5. Here, we assume for simplicity that h\ — h\. 
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Following PJ, we rearrange the singlet KK-Weyl spinors £ an d Xn> sucn that the smallest 
diagonal entry of the KK neutrino mass matrix Ai KK in (I2.16J1 is |e| = min (jM — j^M < 1/ (2R), 
for a given value k = k . In this newly defined basis, the effective kinetic Lagrangian (j2.15j) 
becomes 

(2.17) 



kin 



M 



KK 



where ~§> u is the reordered (4-component) Majorana-spinor vector 



and A4„ the corresponding KK neutrino mass matrix 









(2.18) 



M 



KK 






m 


m 


m 


m 


m 


m 


e 














m 
















m 
















m 
















m 

















V 



(2.19) 



The eigenvalues of A4^ K can be computed from the characteristic eigenvalue equation 
det (A4^ K — Al) = 0, which is analytically given by 

,2 



n ^ 



n=0 



^2 



1 + 



A 



m 



E 



i 



(A - e)- - § 



0. 



(2.20) 



Since it can be shown that A — e = ±n/R is never an exact solution to the characteristic 
equation, only the second factor in (|2.2(J|) can vanish. Employing complex contour integration 
techniques, the summation in the second factor in (J2.2Uj) can be performed exactly, leading to 
an equivalent transcendental equation 



A 



irm R cot 



irR(X-e) 



(2.21) 



As was already discussed in if e = 0, (J2.21j) implies that the mass spectrum consists of 

massive KK Majorana neutrinos degenerate in pairs with opposite CP parities. If e — 1/(2/2), 

the KK mass spectrum contains a massless state, which is predominantly left-handed if mR < 1, 

while the remaining massive KK states form degenerate pairs with opposite CP parities, exactly 

as in the e = case. However, if e ^ 0, 1/(2R), the lepton number gets broken. 2 In this case, 

2 Alternatively, the lepton number may also be broken through the Scherk-Schwarz mechanism, where the 
Scherk-Schwarz rotation angle will induce terms very similar to those depending on e 
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there is no massless state in the spectrum, and the above exact degeneracy among the massive 
Majorana neutrinos becomes only approximate, with a mass splitting of order 2e for each 
would-be (e — > 0) degenerate KK pair. 

We now consider an orbifold model, in which the y = brane is displaced from the orbifold 
fixed points by an amount a. Under certain restrictions in Type I string theory J^Hlll], such 
an operation can be performed respecting the Z2 invariance of the original higher- dimensional 
action. In particular, one can take explicitly account of this last property by considering the 
following replacements in the effective Lagrangian (|2.4|) : 



6<% - 0) 
rj5(y - a) 



21 



+ % + a - 2ttR) 
5(y-a) - 5(y + a - 27rR) 



(2.22) 



with < a < tcR and < y < 2nR. It is obvious that a Z 2 -invariant implementation of brane- 
shifted couplings requires the existence of two branes at least, placed at y = a and y = 2-nR — a. 
In addition, we assume that a is a rational number in units of ttR, i.e. 



ixR 



(2.23) 



where r, q are natural numbers. This last assumption has been introduced for technical reasons. 
It enables us to carry out analytically the infinite summations over KK states (see also our 
discussion below). 

Proceeding as above, the effective KK neutrino mass matrix A4^ K for the orbifold model 
with a shifted brane can be written down in an analogous form 

\ 



M 



KK 



( 


m (0) 


m (!) 


m( 


1) 


m (2) 




2) 


m (0) 


£ 


















m (!) 




















mi- 1 ) 








e — 


1 

R 










m (2) 


















m ("2) 
















e — 


2 
R 



[2.2A) 



where 



rn 



(n) 



V 

V2 



hi cos 



k )a 



R 



ho sin 



k Q )a 



R 



m cos 



na 

7F 



^2.25) 



with m = V\J (h\ + h r ^)j2 and 4>h = tan 1 (/i2/^i) + koa/R. As before, we consider an one- 
generation model with hi = hi and h 2 = h e 2l which renders the analytic determination of the 
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eigenvalue equation tractable. We will relax this assumption in Section 5, when discussing 
the compatibility of this model with neutrino oscillation data. Thus, for our one-generation 
brane-shifted model, the characteristic eigenvalue equation reads 



n 

n=0 



A — £ 



n 

#2 



which is equivalent to 



1 + 



A 



1 



A - e 



X-e 



E 



m 



(n)2 



A — £ - 







E 



m 



(n)2 



A 



S — 75 



(2.26) 



(2.27) 



As opposed to the a = case, complex contour integration techniques are not directly applicable 
in evaluating the infinite sum in (|2.27jl . The preventive reason is that the function mS n \ 
analytically continued to the complex n-plane, is not bounded from above as n — > ±ioo, as it 
had to be, because of its dependence on cos(na/R). However, as has been mentioned above and 
discussed further in Appendix A, this difficulty may be circumvented by assuming that a is a 
rational number in units of irR, as stated in (|2.23|) . Under this technical assumption, we carry 
out in Appendix A the infinite sum in (|2.27|) analytically and derive the eigenvalue equation 
for the simplest class of cases, where a = irR/q with r = 1 and q an integer larger than 1, i.e. 
q > 2. More precisely, we find 3 



A 



irm 2 R 



cos 



— a(X — e) 



cot 



irR(X-e) 



— - sm 
2 



2<p h - 2a(A - e) 



(2.28) 



Observe that unless e = 1/(2/2), a = nR/2 and <ph = 7r/4, the mass spectrum consists of massive 
non-degenerate KK neutrinos. However, it can be shown from ()2.28|) that this tree-level mass 
splitting between a pair of KK Majorana neutrinos is generally small for m( n ) 3> 1/R. In 
particular, this tree-level mass splitting is almost independent of a and subleading so as to play 
any relevant role in our calculations. 

At this stage, it is important to comment on taking the limit a = nR/q — > in (|2.28|) . or 
equivalently q — > oo. This limit is not the eigenvalue equation ()2.21|) which is valid for a = 0, 
because of the presence of the extra non- vanishing term that depends on sin(20h) in (|2.28|) . This 
apparent paradox can be resolved by noticing that the existence of this would-be anomalous 
term is ensured only if the brane-shifting a is much larger than the fundamental quantum 
gravity scale Mp, i.e. a ^> 1/Mp. Since Mp represents a natural ultra-violet cut-off of the 
theory, we expect the onset of new physics above the scale Mp, most likely of stringy nature, 
effectively implying that the KK- Yukawa mass terms m*™** are exponentially suppressed or zero 
for KK-numbers n > MpR. As we will explicitly demonstrate in Section 4 (see our discussion 

3 The so-derived formula generalizes the one presented in |2] to include brane-shifting and arbitrary Yukawa- 
coupling effects. 
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in (|4.18j) ). such a truncation of the KK sum at Mp effectively results in a modification of the 
eigenvalue equation ()2.28j) to 



A = m 2 R 1 7r cos 2 



cj) h - a(X - e) 



cot 



tiR (A - e) 



- Si(2aM F ) sin 



2a(A-e) 



(2.29) 



In the above, Si(x) = Jq dt is the integral-sine function. For any finite value of its argument, 
Si(x) can be expanded as 



+°° (_1 r (2n-l) 

For small x, it is Si(x) ~ x, while Si(x) = 7r/2 for x — > oo. Clearly, as long as a > 1/Mp, the 
eigenvalue equations ()2.28|) and ()2.29|) are almost identical, since Si(2aM F ) = 7r/2 to a very 
good approximation. On the other hand, the limit a — » does now smoothly go over to (j2.21j) . 
as it should be. 

Finally, in addition to the aforementioned tree-level mass splitting, one-loop radiative effects 
may also contribute to further increase the mass difference between two nearly degenerate KK 
Majorana neutrinos, if hi and J12 do not vanish simultaneously. The one-loop generated mass 
splitting, however, is expected to be small [5] of order Tix^m^/ (8ir 2 ) ~ 1CT 2 x (M F /M P ) 2 x 
m( n ) < 10~ 2 x Amy, where mr n ) ~ n/R < M F is the approximate mass of the nth KK pair of 
nearly degenerate Majorana neutrinos, and Am( n ) = m(„+i) — m( n ) ~ l/i? is the mass difference 
between two adjacent KK Majorana pairs. Although such a radiatively-induced mass splitting 
may play a significant role for leptogenesis j^j, its effect on the double beta decay amplitude 
is negligible. Therefore, we neglect radiative effects on the KK mass spectrum throughout the 
paper. 



3 RG evolution of neutrino Yukawa couplings 

The RG evolution of the Yukawa couplings in the standard 4-dimensional scenario involving 
sterile neutrinos has been discussed in [21]. Here, we derive the corresponding RG equations 
for the higher-dimensional case. Since the RG evolution equations for h\ and h l 2 will be similar, 
we concentrate only on the former (= h). In such a higher-dimensional scenario, the presence 
of the KK sterile states alters the RG running. The triangle and self-energy diagrams that 
contribute to the running remain the same as in the SM, except that in the higher dimensional 
context, wherever there are internal £ n lines, there is a multiplicative factor t$ = (fiR) 5 Xs, with 
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Xg = 27r 5 / 2 /8T(8/2). The RG equation for the Yukawa coupling h is given by 
dh 3 



16tt 2 



dlnfi 2 



t s {hh f ) h-h (hlh e ) 1 + hTr ( 3h ] u h u + ?,h\h d + h% + t s h ] h 

h(\gi + \g' 2 ), (3.1) 



where g w and g' are the SU(2)^ and U(l)y gauge-coupling constants, respectively. Note that for 
5 = (1), it is Xs = 1 (2). Also, for 5 = 0, tg = 1, the standard RG equation is reproduced 



We now observe that the four- dimensional Yukawa coupling (h) is suppressed with respect 
to the higher-dimensional coupling (h) by means of the relation: h = (M F /Mp) 5 / n °h. Thus, 
even if we consider h(l/R) ~ 1, the four-dimensional h is suppressed by many orders of mag- 
nitude. From Eq. (|3.1|) . it is also obvious that unless ts is large enough to be comparable with 
(Mp/Mp) 6 ^ 9 , the contributions from the top-quark Yukawa coupling or the gauge couplings 
dominate the running, and hence there is no power-law behaviour at lower energies. 

On the contrary, if we go to a very high energy such that we can ignore h t , then the terms 
multiplying ts dominate. In such a case, ignoring the gauge contribution, we can write 

167r2 TT %Uh?. (3.2) 

dmfi 2 

Integrating Eq. (|3.2|) from the scale /z = R^ 1 to fi, we obtain 

w/r) " w(ti - ~Lk itiR)S - (3 - 3) 

In terms of the Yukawa fine structure constant a(fi) = /i 2 (/z)/(47r) of the original 5-dimensional 
Yukawa coupling (h) and for the simple case § = n g , ()3.3j) takes on the form 

1 1 5X S f n 



a(n) a(l/R) 4tt5 \M f 
Clearly, a(fi) — > oo, for a critical scale 



(3.4) 



/ 4tt5 \ 1/6 
toal = Mf [5X s a(l/R)) (3 ' 5) 

Interestingly enough, ()3.5|) implies that the power-law behaviour sets in not just above the 
compactification scale R~ x , as was naively expected 0, but well above the quantum gravity 
scale Mp. On the other hand, requiring that a(Mp) < 1 in (J3.4)) implies that a(l/R) < 0.55 
for (5=1. This last condition assures that our theory remains perturbative up to the quantum 
gravity scale Mp. From our discussion above, it is obvious that power- law effects on the Yukawa 
neutrino couplings can be safely neglected in our analysis. 
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4 Effective neutrino-mass estimates 



In this section, we calculate the 0^/3/3 observable (m) in orbifold 5-dimensional models. This 
quantity determines the size of the neutrinoless double beta decay amplitude, which is induced 
by jy-boson exchange graphs. To this end, it is important to know the interactions of the W 
bosons to the charged leptons I = e, /x, r and the KK-neutrino mass-eigenstates ri( n ). Adopting 
the conventions of [22], the effective charged current Lagrangian is given by 



CZt = -^ W ~" £ [BivJ^PlVi + E B i,n h,P L n {n) ) + h.c, (4.1) 



=e,^t,r 



where g w is the weak coupling constant, Pl — (1 — 7s)/2 is the left-handed chirality projector, 

and B is an infinite dimensional mixing matrix. The matrix B satisfies the following crucial 
identities: 

+ 0O 

BinB^ + £ B ^ B ln = Siv , (4.2) 

n=— oo 

+oo 

B lvi m n B Vn + E B i,n'rn(n)P>l',n = 0. (4.3) 



Equation (|4.2|) reflects the unitarity properties of the charged lepton weak space, and ()4.3j) holds 
true, as a result of the absence of the Majorana mass terms viVii from the effective Lagrangian 
in the flavour basis. For the models under discussion, the KK neutrino masses mi n ) can be 
determined exactly by the solutions of the corresponding transcendental equations. To a good 
approximation, however, these solutions for large n simplify to 4 

m (n) ~ d + £ ■ ( 4 - 4 ) 

This last expression proves to be a good approximation in our estimates. 
According to (jl.lj) . the Oz//5/3-decay amplitude T^pp is given by [TTj : 

\TJfl) 

= — A^gtfC^) , (4.5) 

where A^gtf = A^gt — A4f is the difference of the nuclear matrix elements for the so-called 
Gamow- Teller and Fermi transitions. Note that this difference of nuclear matrix elements 



For |n| > e and n < 0, the KK mass eigenvalues m( n ) are negative. This corresponds to a neutrino with 
positive physical mass |77i( n \| and negative CP parity. One can always take account of the negative CP parity 
by redefining the mixing matrix elements Bi^ n as -B/.-„ — > iBi ^ n , for n > eR > 0. Although we will allow 
negative neutrino masses in our calculations, wc should stress that both approaches are fully equivalent leading 
to the same analytic results. 
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sensitively depends on the mass of the exchanged KK neutrino in a 0uf3f3 decay. Especially if 
the exchanged KK-neutrino mass mu) is comparable or larger than the characteristic Fermi 
nuclear momentum qp « 100 MeV, the nuclear matrix element .Mgtf decreases as l/m 2 n y 
The general expression for the effective Majorana- neutrino mass (m) in (|4.5jl is given by 

( m ) = T~a 7 V 51 B 2 m {n) \M GT F{m {n) ) - M GT F(m u ) ] . (4.6) 

In the above, the first term describes the genuine higher- dimensional effect of KK-neutrino 
exchanges, while the second term is the standard contribution of the light neutrino u, rewritten 
by virtue of (|4.Hjl . Note that the dependence of the nuclear matrix element Ai G TF on the KK- 
neutrino masses mr n ) has been allocated to (m) in ([4.6)1 . The latter generally leads to predictions 
for (m) that depend on the double beta emitter isotope used in experiment. However, the 
difference in the predictions is too small for the higher-dimensional singlet-neutrino models to 
be able to operate as a smoking gun for different 0i//3/3-decay experiments. 



4.1 Factorization Ansatz for analytic estimates 

To obtain analytic estimates that will help us to gain a better insight into the dynamical 
properties of (|4.6|) . it proves useful to approximate the 0i//?/?-decay amplitude %) V pp in (|4.5jl by 
means of the factorizable Ansatz 

% vPP « (™hA MGTF (m u ) + ^ (m- 1 ) A^gtfK) , (4.7) 

m e m e 

where m p is the proton mass, and -Mgtf(^i/) and ■M-GTF{ m p) are the values of the nuclear 
matrix element A^gtf at m u and m p , respectively. In (|4.7J) . the Ov/3/3 matrix element has 
been written as a sum of two terms. The first term, which is the dominant one, accounts for 
effects coming from KK neutrinos lighter than the characteristic Fermi nuclear momentum qp w 
100 MeV. In this kinematic region, the nuclear matrix element A4 G tf is almost independent of 
the KK neutrino mass m( n ). The second term in (|4.7jl is due to KK neutrinos much heavier than 
qp. This is generically a subdominant contribution to T Qv pp, since M. G TF( m P ) A / 1gtf(^ I /)- 

The quantity (m)sA is an approximation of the effective Majorana-neutrino mass (m), which 
is obtained by approximating the nuclear matrix elements M. G r FF{ m {n)) entering (m) in (14.6)1 
by a step function at |m( n )| = g^: 

( v /-MgtfW, for |m (n) | < g F , 
Mgtf(W(„)J - < (4.8) 
1 , for |m (n )| > g F • 
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In what follows, we refer to such an approach to the nuclear matrix elements as the Step 
Approximation (SA). The effective neutrino mass in the SA reads: 

l(q F -e)R] 

(m) S A = B 2 ev m v + B e,n m (n) 

n=-[(q F +e)R] 
+00 +00 

= - E B ln m (n) ~ E B e,-n m (-n) , (4.9) 

n=[(<lF—e)R] n=[(q F +e)R] 

where we used (14.3)1 to arrive at the last equality for the effective neutrino mass. Notice that 
(m) is not zero, simply because the sum over the KK neutrino states is truncated to those with 
a mass |wi( n )|, |j7i(- n )| < Qf- 

Correspondingly, the effects of the heavier KK neutrinos, with masses mr n ) ~ Qf, have been 
taken into account in the factorizable Ansatz ()4.7jl by means of the inverse effective neutrino 
mass (m" 1 ). This newly introduced quantity is given by 

+00 +00 

(™~ 1 > = E B ln^n) + E B l-U^\-n) ■ (4-10) 

n=[(q F -e)R] n=[(q F +e)R] 

The factorizable form (|4.5jl of the matrix element constitutes a good approximation except for 
the isolated region where |m,( n )| ~ qF ~ 100 MeV. Nevertheless, the effect of the KK neutrinos 
on the effective neutrino mass is cumulative [H] due to a sum of an infinite number of states, 
since each KK state has either a tiny Majorana mass or a very suppressed mixing with the 
electron neutrinos. Therefore, we expect that excluding this isolated region of KK-neutrino 
contributions around qp will not alter quantitatively our results in a relevant way. 

We will now rely on (J4.9)) to estimate the effective neutrino mass (m)sA m different settings 
of 5-dimensional orbifold models discussed in Section 2. To begin with, let us consider a simple 
orbifold model, with e 7^ and e 7^ 1/(2R). In addition, we consider the case a = 0, namely 
we take the brane to be located at the one of the two orbifold fixed points. Like the neutrino 
masses, the mixing- matrix elements B eu and B Sjn can also be computed exactly [2]: 

B ev = 1, B e>n = (4.11) 
where the squares of the normalization factors M and N( n ) are given by 

^ = 1 + £ t ^— tt , K> = 1 + E t ^— ^ ■ (4.i2) 

n=-oo ^£ - m y + ^ 1 fc=-oo ( £ - m (n ) + 

Applying complex integration methods for convergent infinite sums, the squared normalization 
factor J\f 2 can be calculated to give 

N 2 = I + . 2 f^ = 1 + n 2 m 2 R 2 + ^ • (4-13) 
sin [TrH(Tn u — e)\ m z 
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In obtaining the last equality in ([4.13)1 . we used the eigenvalue equation ([2.21)1 for A = m v . 
From ([4.11)1 and ()4.13|) . we immediately see that if mR -C 1 and m u <C m, it is B eu « 1 
and hence the lightest neutrino state is predominantly left-handed. For the calculation of the 
effective neutrino mass, we need 



,. 9 ir 2 m 2 R 2 9 9 9 m 2 n) (^ + e ) 2 

M? n) = 1 + , 2r = 1 + 7T 2 m 2 R 2 + « Vfi 9 ^ , (4.14) 

v ; sm [TrH[mt n ) — £)\ m m 



where the last approximate equality in (|4.14|l corresponds to a large n. In Appendix B, we 
show that the KK neutrino masses derived from (|2.21j) and the mixing-matrix elements given 
in (|4.11|) satisfy the sum rules given by the identities ()4.2|) and ([4.3)1 . 

Based on (|4.9j) . we will now perform an estimate of the effective neutrino mass in the simple 
orbifold model mentioned above. Plugging the value of B en = 1/Ni n ) into ([4.9)1 . we may 
estimate the effective neutrino mass in the SA through the following steps: 

00 / 1 1 \ 

( m )sA = E ~ I n + " n + O 





R ' 

1 A _ ^ d1 J9f- £ 



(4.15) 



In arriving at the last equality in ([4.15)1 . we approximated the sum over the KK states by an 
integral, and used the fact that e/qp <C 1. Since 2eR < 1, we can estimate that for m — 10 eV, 
( m )sA ~ 10 -6 eV, which is undetectably small. 

The above large suppression of the effective neutrino mass (m)sA is a consequence of the 
very drastic cancellations due to KK neutrinos with opposite CP-parities. However, we might 
be able to overcome this difficulty by arranging the opposite CP-parity KK neutrinos to couple 
to the electron and W boson with unequal strength. In fact, this is what happens in orbifold 
models automatically, if the y = brane is shifted to y = a 7^ 0. In this case, the mixing-matrix 
elements B eu and B e>n are given by the inverse of M and Af( n ) respectively; but now for the 
shifted brane, A/"( n ) is given by 

+°° cos 2 (*2 - 6 h ) f»+ £ )2 

"»- 1 + *& (.-^!ly ~*#T^- (416) 

where the second approximate equality in ([4.16)1 corresponds to large n. 
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By analogy to (J4.15|) . we may compute the effective Majorana-neutrino mass for the the 
brane-shifted scenario (a ^ 0) as follows: 



/\ 2n f , r os2 (f zjjO cos2 (f +< n n 

(m)sA W - mjR 7 H n + eR ~ n- eR ) 

q F R 

MpR / 9 \ 

. , , „ f an /2na\ / em R \ . , . 

= -sm(2« m ^ / -s,n — + O — . (4.17) 

In the second step, we have truncated the upper limit of the integral at the fundamental quan- 
tum gravity scale Mp. The scale Mp represents a natural ultra-violet cut-off of the problem, 
beyond of which the onset of string-threshold effects are expected to occur. The last result 
in 1)4.17)1 can now be expressed in terms of the integral-sine function Si(x) = Jq dt^f^. Thus, 
the effective neutrino mass can be given by 



(ra)sA ~ - sm(2(j) h ) m 2 R 



Si(2aM F ) - Si(2agi 



+ (2£5). ,4,8, 



Notice that for a fixed given value of Mp, the analytic expression ()4.18j) for the effective neutrino 
mass goes smoothly to ()4.15|) in the limit a —>■ 0, as it should be. In order that the prediction 
for neutrinoless double beta decay effects is at the level reported recently [T3], we only need to 
have: <ph ~ ±vr/4 and 1/Mp <d< l/(2qp), i.e. the brane is slightly displaced from its origin. 
For instance, if a ~ 1/(3^), m = 10 eV and 1/R = 300 eV, we find that (m)sA is exactly at 
the observable level, i.e. (m)sA ~ 0.4 eV. 

It is now interesting to give an estimate of the inverse effective neutrino mass (m _1 ) in 
the orbifold model with a shifted brane (a 7^ 0). The quantity (m _1 ) can be approximately 
calculated as follows: 

+ r f COS 2 (?§ - (f) h ) COS 2 f + <P h ) \ 

(m- 1 ) « m 2 R 3 / dn ( . U J" J - -. U J?' ) (4.19) 



/ , n 9^-? f dn f2na\ 3 . . , . f dn f2na\ 
sin(20 h ) m 2 R 3 J — sin f_j + - cos(20 /l ) m 2 ^ 4 y — sin ( — J 



The RHS of the last equality in (|4.19|) can be written down in a lengthy expression in terms of 
the integral-sine, integral-cosine and known trigonometric functions. For example, for (f>h = 7r/4, 
(m _1 ) is given by 



(m x ) m 2m 2 R 



Si (2aqp) — — ) — — 5- sin(2agi?) — cos(2aqp) 

2 J Aqt 2q F 



em 2 R 



2q 3 



F 

(4.20) 
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For the specific model considered above, with m = 10 eV, 1/R = 300 eV and a = l/(3q F ), 
we find that (m _1 ) < 10~ 5 TeV -1 . Hence, the above exercise shows that the contribution from 
(jrT 1 ) to the double beta decay amplitude ()4.5|) is subdominant; it gets even more suppressed 
for a <C 1/qF- 

4.2 Numerical evaluation 

To obtain realistic predictions for the double beta decay observable (m), one has to take into 
account the dependence of A^gtf on the KK neutrino masses mi n )- To properly implement 
this m( n )-dependence in our extractions of the effective Majorana mass (m) from the different 
nuclei, we have used the general formula (|4.6|) . where the infinite sum over n has been truncated 
at l^-maxl = MpR, namely at the quantum gravity scale Mp. Notice that the general formula 
for (m) in ()4.6|) includes the contributions from the KK neutrinos heavier than qp, described 
by the inverse effective neutrino mass (m- 1 ) in (J4~20j) . 

In Table ^ we present numerical values for the difference of the nuclear matrix elements, 
-Mgtf = .Mgt _ M-Fi as a function of the KK neutrino mass mi n )- Our estimates are obtained 
within the so-called Quasi-particle Random Phase Approximation (QRPA) [2UEE]- Here, we 
should note that the numerical values for the nuclear matrix element of 100 Mo exhibit some 
instability due to its sensitive dependence on the particle-particle coupling g PP within the 
context of the QRPA. In addition, we should remark that in our numerical evaluation of (m), 
the nuclear matrix elements A^gtf have been interpolated between the values given in Tabled 

In Tabled we show numerical values for the effective Major ana-neutrino mass (m) as derived 
for different nuclei in a 5-dimensional brane-shifted model, with m = 10 eV, 1/R = 300 eV, 
e = 1/(4R), <ph = — 7r/4 and Mp = 1 TeV. In addition, we have varied discretely the brane- 
shifting scale 1/a from 0.05 GeV up to values much larger than Mp. The first column in TableEl 
give the predictions obtained in the SA for the nuclear matrix elements. The SA is closely 
related to our approximative method followed above, leading to results that are in a very good 
agreement with 1)4.18)1 . Remarkably enough, even the change of sign of (m)sA at 1/a ~ 0.1 GeV 
in Table El can be determined sufficiently accurately by analyzing the multiplicative expression 
7r/2 — Si(2aqp) in ()4.18j) . which oscillates around n/2 [26J, for 1/a < 0.1 GeV. Analogous remarks 
can be made for the inverse effective neutrino mass (m" 1 ) in (|4.20p . 

As can be seen from Table El the deviation between the SA and the one based on the 
general formula ()4.6)) is rather significant if a is close to 1/qF due to the non-trivial nuclear 
matrix element effects mentioned above and due to heavier KK-neutrino effects coming from 
(rrT 1 ). However, for smaller values of a, i.e. for a < l/(3qp), the agreement between the 
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m (n) [ lvl c V J 






M 


GTF(ra(n)) 










76 Ge 


82 Se 


100 Mo 


116 Cd 


< 1 


4.33 




4.03 




4.86 




3.29 




10 


4.34 




4.04 




4.81 




3.29 




10 2 


Q HQ 




Z.oZ 




3.31 




z.lo 




10 3 


1.40 x 10" 


-l 


1.25 x 10- 


-1 


1.60 x 10" 


-l 


9.34 x 10" 


-2 


i n4 
1U 


1.39 x 10" 


-3 


1.24 x 10- 


-3 


1.60 x 10" 


-3 


9.26 x 10" 


-4 


10 5 


1.39 x 10" 


-5 


1.24 x 10- 


-5 


1.60 x 10" 


-5 


9.26 x 10" 


-6 


1U 


1.39 x 10- 


-7 


1.24 x 10- 


-7 


1.60 x 10" 


-7 


9.26 x 10" 


-8 


10 7 


1.39 x 10- 


-9 


1.24 x 10- 


-9 


1.60 x 10" 


-9 


9.26 x 10- 


10 


m (n) [ 1V1( = V J 






M 


GTF(^(n)) 










128 Te 


130 Te 


136 Xe 


150 Nd 


< 1 


4.50 




3.89 




1.83 




5.30 




10 


4.52 




3.91 




1.88 




5.45 




10 2 


3.19 




2.79 




1.48 




4.24 




10 3 


1.46 x 10- 


-1 


1.29 x 10- 


-1 


7.07 x 10" 


-2 


2.02 x 10" 


-1 


10 


1.46 x 10- 


-3 


1.28 x 10- 


-3 


7.04 x 10" 


-4 


2.02 x 10" 


-3 


10 5 


1.46 x 10- 


-5 


1.28 x 10- 


-5 


7.05 x 10" 


-6 


2.02 x 10" 


-5 


10 6 


1.46 x 10- 


-7 


1.28 x 10- 


-7 


7.05 x 10" 


-8 


2.02 x 10" 


-7 


10 7 


1.46 x 10- 


-9 


1.28 x 10- 


-9 


7.05 x 10- 


10 


2.02 x 10" 


-9 



Table 1: QRPA estimates of the relevant combination of nuclear matrix elements, A^gtf = 
A^gt — M-f, as o, function of the KK neutrino mass rri( n ). 

effective neutrino mass computed in the SA and the general formula ()4.6|) is fairly good. In 
this kinematic regime, the inverse effective neutrino mass (m -1 ) becomes rather suppressed 
according to our discussion in ()4.20|) . Our numerical estimates in the last column of Table El 
offer firm support of this last observation. Thus, the main contribution to (m) originates from 
KK neutrinos much lighter than qp. Consequently, within the 5-dimensional brane-shifted 
model, we have numerically established a sizeable value for (m) in the presently explorable 
range 0.05-0.84 eV. Finally, for very small values of a, i.e. for a -C 1/Mp, we recover the 
undetectably small result ()4.15|) for the unshifted brane a = 0. 
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l/a 


(m) [eV] 


/ 1 \ 

{m ) 


[GeV] 


SA 


76 Ge 


82 Se 


100 Mo 


116 Cd 


128 Te 


130 Te 


136 Xe 


i50 Nd 


[TeV^ 1 ] 


0.05 


-0.062 


0.009 


0.010 


0.016 


0.012 


0.009 


0.008 


-0.004 


-0.004 


6.2xl0~ 6 


0.1 


-0.012 


0.052 


0.054 


0.061 


0.062 


0.052 


0.050 


0.025 


0.026 


-3.6 xl0~ 6 


0.2 


0.208 


0.096 


0.100 


0.109 


0.114 


0.097 


0.094 


0.058 


0.061 


-1.3xl0~ 5 


0.3 


0.307 


0.123 


0.128 


0.136 


0.143 


0.124 


0.121 


0.082 


0.086 


-1.2xl0~ 5 


1 


0.457 


0.271 


0.275 


0.280 


0.287 


0.272 


0.269 


0.241 


0.243 


-5.7xl0~ 6 


10 


0.516 


0.493 


0.493 


0.494 


0.495 


0.493 


0.493 


0.489 


0.489 


-6.6 xHT 7 


10 2 


0.515 


0.513 


-6.7xl0" 8 


10 3 










0.535 










-6.7xl0~ 8 


10 4 










0.066 










-6.9xl0~ 10 


10 10 










< io- 6 










0. 



Table 2: Numerical estimates of (m) for different nuclei in a 5- dimensional brane- shifted model, 
with m = 10 eV, 1/R = 300 eV, e = 1/(47?), <f> h = -vr/4 and M F = 1 TeV. The first column 
exhibits the numerical values for (m) in the Step Approximation ( SA ) for the nuclear matrix 
elements, while the last column shows the results for the inverse effective neutrino mass (m _1 ). 



4.3 (m) and the neutrino mass scale 



Apart from explaining the recent excess in 0u(3(3 decays, the 5- dimensional model with a 
small but non-vanishing shifted brane exhibits another very important property. The effective 
Ma jorana- neutrino mass (m) can be several orders of magnitude larger than the light neutrino 
mass m u , for certain choices of the parameters e and <ph- To understand this phenomenon, let 
us first consider the eigenvalue equation ()2.27|) for A = m u , written in the form: 

00 m (")2 

Notice that (j4.21|) constitutes an excellent and very practical approximation of the neutrino- 
mass-mixing sum rule, when the small m^-dependence in the infinite sum over the KK neutrino 
states is neglected and the approximate formulae ()4.4)1 and 1)4.16)1 for the KK masses mr n ) and 
mixing-matrix elements B e<n , along with B ev = 1, are substituted in (|4.3|) . Then, the infinite 
sum over KK neutrino states can be performed with the help of (J2.28)) for rational values of 
a in irR units. Especially for a = irR/q with q being an integer much larger than 1, i.e. for 
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1/M F <a< 1/<?f, the light neutrino mass m u is given by 



m,, « — nm 2 R 



cos 2 ft cot (vr .Re) + -sin(20 h )]. (4.22) 

It is now easy to see that the light neutrino mass m v can be very suppressed for specific values 
of 4>h and e. For instance, one obvious choice would be <ph ~ — vr/4 and e ~ 1/(4R). On the 
other hand, the effective neutrino mass (m)sA is determined by the second sine-dependent term 
in (|4.22j) (cf. (|4.18jl ). which is induced by brane- shifting effects. Unlike the suppressed light 
neutrino mass m u , the effective neutrino mass (m) can be sizeable in the observable range 0.05- 
0.84 eV. This loss of correlation between the quantities (m) and m v is a rather unique feature of 
our higher- dimensional brane-shifted scenario. As we will discuss in the next section, the above 
de-correlation property plays a key role in our model-building of 5-dimensional brane-shifted 
scenarios that could explain the neutrino oscillation data. 



5 Atmospheric and solar neutrino data 

Atmospheric and solar neutrino data [2*7|I28(I29( , together with information from laboratory ex- 
periments, such as the CHOOZ experiment [3U], are very crucial for a given higher- dimensional 
singlet-neutrino model to qualify as viable. In particular, the latest SNO results [21] appear 
to disfavour large components of sterile neutrinos, indicating a preference among the different 
solutions to the solar and atmospheric neutrino puzzles for those involving transitions between 
almost active neutrinos. 5 To account for this experimental indication, we assume that the 
compactification scale 1/R and the lepton-number- violating bulk parameter e are much larger 
than the KK Dirac mass terms in (12.25(1 . 

In the following, we shall explicitly demonstrate that our 5-dimensional brane-shifted model 
with only one bulk neutrino is able to fully explain the neutrino oscillation data. Specifically, 
we will show that the preferred solar Large Mixing Angle (LMA) and atmospheric solutions, 
which both require large v^-v^ and u^-i/j- mixings, can be realized within our 5-dimensional 
model. These particular solutions are allowed, only if the differences of the squares of the light 
neutrino masses lie in the ranges: 

1.8 xl0~ 3 < Am 2 tm [eV 2 ] < 4.0 x 10~ 3 , 2.0 x 10~ 5 < Am 2 Q [eV 2 ] < 2.0 x 10~ 4 , (5.1) 

with Am 2 tm = m 2 3 — m 2 2 and AtUq = m 2 2 — vr? vx . According to the usual conventions, the 
physical light neutrino masses m Ul , m V2 and labelled in increasing hierarchical order, 

i.e. m ui < m U2 < m U3 . 

5 A recent study [HI] seems to suggest that the active neutrino component in the solar neutrinos has to be 
larger than 86% at 1 a CL. A loophole may exist for atmospheric neutrinos, see |82j . 
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To start with, let us consider the weak basis in which the charged lepton mass matrix is 
diagonal. Then, in the three-generation brane-shifted model, the KK-Dirac Yukawa terms are 
given by the 3-vectors 



m 



(n) 



m cos 



cos 



where 



m 



v 

V2 



-m V(M) 2 + (h l 2 . 



na 
R 


4>e)\ 


na 




R 


<t>„) 


na 


<t>r)j 


R 


h = 


tan - 



k a 
~R 



(5.2) 



(5.3) 



with / = e,/i, r. Given our assumption that e,l/R 3> m l , the KK neutrinos can now be 
integrated out. Analogously with ()2.27|) . the effective light neutrino mass matrix Ai u can be 
computed by 

+00 (n) m (n)T 

M v = - £ ■ (5-4) 



R 



+ e 



Following the same line of steps as in Appendix A, one is able to analytically carry out the 
infinite sum in (|5.4J) for the phenomenologically interesting case of a = nR/q, with q being an 
integer much larger than 1. In this limit, we obtain the novel trigonometric mass texture: 



Mr,, 



TxRm l m l 



cos <pi cos (pi cot^TcRe) + - sin(0; + <pi. 



(5.5) 



with I, I' = e, /i, r. The effective neutrino mass matrix ()5.5|) consists of two terms: (i) the cosine- 
dependent term that arises from the lepton-number-violating bulk mass M (or equivalently 
e) and (ii) the sine-dependent term which is due to lepton-number violation in the effective 
Yukawa couplings and is caused by slightly shifting the brane from the orbifold fixed points. 
The occurrence of the second brane-shifting mass term is always ensured as long as a > 
1/Mp. Without the presence of this brane-shifting- induced term, the effective neutrino mass 
matrix ()5.5|) is of rank 1, leading to two massless neutrinos. This last fact is very undesirable, 
as it would be very difficult to explain both solar and atmospheric neutrino data with only 
one non-trivial difference of neutrino masses in the frequently-discussed scenario without brane 
shifting. 

As has been discussed in Section 4, however, even a small amount of brane shifting may in- 
duce sizeable lepton-number-violating Yukawa interactions. The latter generate brane-shifting 
mass terms that break the rank-1 structure of the effective neutrino mass matrix Ai v '. The 
resulting Ai u in (15.5(1 exhibits a novel trigonometric structure that can predict hierarchical 
neutrinos with large v v ,—v T and v^— v e mixings to explain the atmospheric and solar neutrino 
anomalies, along with a small v e -v T mixing as required by the CHOOZ experiment [HD]- At this 
point, it is important to stress that the effective neutrino mass (m) entering the O^/3/3-decay 
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amplitude gets fully decoupled from the neutrino-mass matrix element M. v ee - According to our 
discussions in Section 4 (cf. ()4.18J1 ). the effective neutrino mass for the three-generation case is 
given by 

a R * AC- (5-6) 



(m) « -- sin(20 e 



7r (m 



It is important to recall again that unlike A4^ e , KK neutrinos heavier than the Fermi nuclear 
momentum qp do not contribute significantly to (m), leading to the loss of correlation between 
(m) and M. u ee . The latter is a distinctive feature of the KK-neutrino dynamics. This de- 
corellation between (m) and M. u ee permit us to consider the interesting case |(m)| ^> l-Mf^l, 
for all I, I' = e, /x, r. Such a realization enables us to accommodate a sizeable positive signal of 
Ou/3/3 decays together with the present neutrino oscillation data. 



To realize the aforementioned hierarchy \{m)\ 3> |.Mf^|, we assume that all phases <pi are 
close to — tt/4. For concreteness, we adopt the following scheme of phases: 

<f>l = - j + S t , ttRs = j - 5 e . (5.7) 

where Si, S E <C 1. Our choice of phases has been motivated by the fact that the above- described 
de-correlation between (m) and Ad^e becomes fully operative in this case. To implement the 
CHOOZ constraint in our model-building, we require that M u eT = -M^e = 0- This last constraint 
implies that 

25 £ = -5 e -5 T . (5.8) 

Moreover, without loss of generality within our phase scheme, we may take 5^ = 0. Under these 
assumptions, the light neutrino-mass matrix takes on the simple form 



M " = if 



\ 



/ m e2 (5 T — 5 e ) m e m fM 5 T 

m fl m e 6 T m fl 2 (5 e + S T ) m fl m T 5, 
\ m T m tM 5 e m r2 (5 e — 5 T ) J 



(5.9) 



Let us now consider the following numerical example: 



Sr 



mr 



— « 1.468 
This leads to the neutrino mass matrix: 



7rm e2 R 



= 25, 
m T 

— « 2.542 



( -1 1.47 \ 

1.47 6.46 7.46 
\ 7.46 6.46 j 



(5.10) 



(5.11) 



Notice that all elements of the neutrino-mass matrix Ai u in (|5.11|) can be suppressed by choosing 
a small value for the factorizable parameter 5. In our numerical example, the neutrino mass 
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matrix (|5.11|) can be diagonalized through v )l -v T and v^-Vp mixing angles close to tt/4, whereas 
the v e -v T mixing angle is small, below 0.1. In addition, its mass-eigenvalues are approximately 
given by 

(AHdiag « 57Tm e2 R{0, 1,7). (5.12) 

Assuming that m e = 10 eV and 1/R = 300 eV for a successful interpretation of the recent excess 
in 0uf3f3 decays, then it should be 5 = (6-9) x 10~ 3 to accommodate the neutrino oscillation 
data through the LMA solution. In particular, we obtain the neutrino-mass differences: 

Am atm ~ (2 - 4) x 10~ 3 eV 2 , Ami ~ (4 - 8) x 10~ 5 eV 2 (5.13) 

These results are fully compatible with the currently preferred atmospheric and solar LMA 
solutions to the neutrino anomalies. 

In our demonstrative analysis carried out in this section, we have not attempted to fit 
the results of the Liquid Scintillator Neutrino Detector (LSND) as well jHHl- In principle, our 
brane-shifted 5-dimensional models are capable of accommodating the LSND results through 
active-sterile neutrino transitions. In this case, however, the lowest-lying KK singlet neutrinos 
should be relatively light. As a result, they cannot be integrated out from the light neutrino 
spectrum, thereby leading to a much more involved effective neutrino-mass matrix. A complete 
study of this issue, including possible constraints from the cooling of supernova SN1987A [BUM], 
is beyond the scope of the present paper and may be given elsewhere. 



6 Conclusions 



We have studied the model-building constraints derived from the requirement that KK singlet 
neutrinos in theories with large extra dimensions can give rise to a sizeable 0^/3/?-decay signal 
to the level of 0.4 eV reported recently. Our analysis has been focused on 5-dimensional 
S 1 1 'Z-i orbifold models with one sterile (singlet) neutrino in the bulk, while the SM fields are 
considered to be localized on a 3-brane. In our model-building, we have also allowed the 3- 
brane to be displaced from the S 1 /Z 2 orbifold fixed points. Within this minimal 5-dimensional 
brane-shifted framework, lepton-number violation can be introduced through Majorana-like 
bilinears, which may or may not arise from the Scherk-Schwarz mechanism, and through lepton- 
number-violating Yukawa couplings. However, lepton-number-violating Yukawa couplings can 
be admitted in the theory, only if the 3-brane is shifted from the S l /Z 2 orbifold fixed points. 
Apart from a possible stringy origin [20], brane-shifting might also be regarded as an effective 
result owing to a non-trivial 5-dimensional profile of the Higgs particle [HE] and/or other SM 
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fields |3b| l3*7] that live in different locations of a 3-brane with non-zero thickness which is 
centered at one of the S 1 /Z 2 orbifold fixed points. 

One major difficulty of the higher- dimensional theories is their generic prediction of a KK 
neutrino spectrum of approximately degenerate states with opposite CP parities that lead to 
exceedingly suppressed values for the effective Majorana- neutrino mass (m). Nevertheless, 
we have shown that within the 5-dimensional brane-shifted framework, the KK neutrinos can 
couple to the W ± bosons with unequal strength, thus avoiding the disastrous CP-parity can- 
cellations in the O^/3/3-decay amplitude. In particular, the brane-shifting parameter a can be 
determined from the requirement that the effective Majorana mass (m) is in the observable 
range jT3j: 0.05-0.84 eV. In this way, we have found that 1/a has to be larger than the typical 
Fermi nuclear momentum qp = 100 MeV and much smaller than the quantum gravity scale 
Mp, or equivalently 1/Mp <C a < 1/qF- 

An important prediction of our 5-dimensional brane-shifted model is that the effective 
Majorana- neutrino mass (m) and the scale of light neutrino masses can be completely de- 
correlated for certain natural choices of the Majorana-like bilinear term e and the original 
5-dimensional Yukawa couplings h\ and h l 2 in ()2.4|) . For example, if e ~ 1/(AR) and h\ ~ —h l 2 , 
we obtain light-neutrino masses that can be several orders of magnitude smaller than (m). 
Nevertheless, it is worth mentioning that if future data did not substantiate the presently re- 
ported 0vj3P excess, the above model-building conditions would then need be modified. Such 
a possible modification would not jeopardize though the viability of our brane-shifted scenario. 
Indeed, if the upper limit on the effective neutrino mass became even lower and lower, this 
would imply that the above decorrelation property is less and less necessary. 

Another important prediction of the 5-dimensional brane-shifted model with only one bulk 
sterile neutrino is that the emerging effective light-neutrino mass matrix does no longer possess 
the rank-1 form, as opposed to the brane-unshifted a = case. As we have shown in Section 5, 
the above properties of the brane-shifted models are sufficient to explain, even with only one 
neutrino in the bulk, the present solar and atmospheric neutrino data by means of oscillations 
of hierarchical neutrinos with large Vg-v^ and maximal v p -v T mixings. In particular, neutrino- 
mass textures can be constructed that utilize the currently preferred LMA solution, where the 
v e -v T mixing is small in agreement with the CHOOZ experiment. 

Although a sizeable 0z//3/3-decay signal can be predicted within our brane-shifted 5- 
dimensional models, the above-described de-correlation property between (m) and the actual 
light neutrino masses suggests, however, that it is rather unlikely that such a signal be accompa- 
nied by a corresponding signal in Tritium beta-decay experiments. For example, the KATRIN 
project [38J has a sensitivity to active neutrino masses larger than 0.35 eV at 95% CL, and so 
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it can only probe the existence of light neutrinos much heavier than those considered in our 
5-dimensional models. Finally, the brane-shifted models under study also have the potential to 
accommodate the LSND results by virtue of active-sterile neutrino oscillations. In this case, 
the lowest-lying KK-neutrino states will contribute to the effective light neutrino-mass matrix, 
giving rise to more involved mass textures. In this context, it would be very interesting to 
investigate the question whether a simple higher- dimensional model accounting for all the ob- 
served neutrino anomalies can be established. We plan to address this interesting question in 
the near future. 
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A Eigenvalue equation 



Starting from (J2.27)) . we will derive here the transcendental eigenvalue equation ()2.28j) . for the 
simplest class of brane-shiftings with a = nR/q, where r = 1 and q is an integer larger than 1, 
i.e. q > 2. Then, the eigenvalue equation ()2.27|) can be equivalently written as 

<?-l oo rn (qk+l)2 q-l oo i 

A = E E ™ -<*±l = S - ( ° 2 E A _/_^ > (a- 1 ) 

1=0 k=-oo A E R 1=0 k=-oo A E R 

where we have used the periodicity property (m"' the second step of (|A.lj) . 

In fact, it is this last periodicity property of the KK- Yukawa terms that we wish to exploit 
here to carry out analytically the infinite sums in (|A.1|) . which has forced us to introduce the 
technical constraint (|2.2Hjl . namely that a/(irR) is a rational number. Now, the individual l- 
dependent infinite sums over k in (|A.1|) can be performed independently, using complex contour 
integration techniques. In this way, we obtain 



A = - nm 2 R ^2 c °s 2 ( 4>h — — ) cot 



1 DM \ l?l 

- ixR (A — e) — — 

q q 



(A.2) 



Our next task is to carry out the summation over I in (|A.2j) . For this purpose, we express 
the RHS of (|A.2|) entirely in terms of sine and cosine functions by factoring out the common 
divisor, i.e. 

A = — x £ cos 2 - - cosf - - - ) II ^{ ~ ~ — )'( A - 3 ) 

/nsin(^) «=o V V ^ « J {Z% V ? ? J 

with 6 = 7i R (A— e). To further evaluate (jA.3|) . we exploit the following trigonometric identities: 6 
> cos IT sm = - '—— q cos 6 , (A.5) 



^ / 2/tt\ /0 Ztt\ 9 ~ x (Q mix 
> cos 20 ft cos . I sin 



-Di- 1 / o — 2 

g cos 20 h + 9 ) . (A.6) 



2 g-l - V ™ 5 



6 The proof of these identities is rather lengthy and relies on the particular properties of the g-roots of the 
unity, i.e. the roots of the equation z q — 1. Specifically, we used the basic property of the unit roots that their 
sum and the sum of their products are zero, while their total product is (— l) 9-1 . 
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With the help of (lA~4l) - (|A~5j) . we arrive at the transcendental eigenvalue equation 

nm 2 R( r /x , cos[20 h + ^7ri2(A-e)] 1 

A = —-{ cot 7ri?(A-£) + L , r * 1 • (A.7) 

z { sin 7r_R (A — e) J 

If we replace q with nR/a in (jA.7|) . we arrive after simple trigonometric algebra at the tran- 
scendental eigenvalue equation ()2.28|) . Although we focused our attention on the simplest class 
with a = irR/q, we should remark that our methodology described above can apply equally 
well to the most general case where the brane-shifting a is any rational number r/q in nR units. 



B Sum rules 

In this appendix, we will show that the KK-neutrino masses determined by the roots of (j2.21|) 
and the mixing- matrix elements given in (j4.11|) satisfy the sum rules f!4.2|) and ()4.3|) . For 
simplicity, we consider the case a = 0. However, our considerations carry over very analogously 
to the case a = irR/q ^ 0, where q is an integer larger than 1. 

Let us first consider (j4.2j) for I = V = e. We will then prove that 

N 

\B ev \ 2 + lim Yl \ B e,n? = 1 • (B.l) 

N ^°° n=-N 



Our proof will rely on Cauchy's integral theorem. Thus, the LHS of (|B.lj) can be expressed in 
terms of a complex integral as follows: 

N x , 1 N 

\B ev \ 2 + lira £ \B e ^ = — lira A dz _ — + £ 



a^oo n ±^ N ' 2-ni n^oo J Cn \z - m u n ±^ N z - m (n) 

1 

x 



1 + ir 2 m 2 R 2 / sin 2 [irR(z-e)} 

If 1 

= lim i dz f . (B.2) 

2iri n-^oo Jc N z — irm 2 R cot[n R(z — e)\ 

In deriving the second equality in ()B.2|) . we have noticed that for z in the vicinity of the pole, 
e.g. for z w m( n ), it is 

f Ti 2 m 2 R 2 1 

z - nm 2 R cot[rcR(z - e)} « (z - m (n) ) \ 1 + . 2r -r > . (B.3) 

I sin [kR{z — e)\ J 

Such a substitution is only valid under complex integration, provided there are no singularities 
of the complex function cot[nR(z — e)\ on the contour Cjv- For this purpose, we choose our 
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contours to be circles represented in the complex plane as 

ZN = (JLt^L + e . (B.4) 

Then, it can be shown that on the complex contours z = zn, | cot ttR(zn — e)| is bounded 
from above by a constant independent of N. Thus, on Cat the last integral in (jB.2|) may be 
successively computed as 

lim </> dz 



2ni n^oo Jc N z — nm 2 R cot[ir R(z — e)] 



1 lim r d e t{zN ~ £) 



2ni n~*oc Jo z N - nm 2 R cot[7r(iV + |) e ie ] 



1 in 7cm 2 R cot[vr(iV + ±) e ie ] - e 
— lim / dO ; — 7 rr — — 

2tt n^oo Jo Zn - irm 2 R cot[7r(iV + |) e %e \ 



The second term in the last equality of ()B.5|) vanishes in the limit N — > oo or equivalently when 
zn is taken to infinity in a discrete manner as prescribed by (|B.4j) . Thus, the complex integral 
in the last equality of (|B.2)1 is exactly 1, which proves the unitarity sum rule (IB. II) . 



In the remainder of the appendix, we will prove the neutrino-mass-mixing sum rule: 

N 

N—>oo 



B 2 eu m u + lim ]T B 2 n m {n) = 0. (B.6) 



n=-N 



In our proof, we will follow a path very analogous to the one outlined above for showing (|B.1|) . 
Thus, the LHS of ()B.6j) may be expressed in terms of a complex integral as follows: 



N If z 

B 2 eu m v + lim^ Bl n m (n) = — Km ^ dz - —^—^-^ - . (B.7) 



n=-N 



2iri n^oo Jc N z — irm 2 R cot [irR(z — e)] 



Evaluating the complex integral on the contours Cn defined by flB.4|) yields 
lim </" dz 



2iri n^oo Jc N z — irm 2 R cot[ir R(z — e)] 



/•a* i( ZN - £ ) nrn ?R C ot[vr(iV + ±) e ie ] 
lim d6 - 



2ni n->oo Jo z N - irm 2 R cot[vr(iV + \)e id ] 

m 2 R Jiiri^ | dO cot[vr(iV + §) e ie ] + 0{l/z N ) \ . (B.8) 



l n r r 2 * 

2 

Similar to the second term in the last equality of (|B.8|) . which goes to zero for N — > oo, the 
first term vanishes as well after integration over 0. This can be readily seen by exploiting 



29 



respectively the periodic and antisymmetric properties of the integrand with respect to 6 and 
its argument: 

dQ cot[7r(iV+ \)e ie \ = I* d6 cot[n{N + \) e ie ] + d6 cot[K{N + \) e ie ] 

JO J J 7T 

d9 cot[n(N +l)e ie ] + F dB cot[vr(iV + |) e l{e+ ^} 



II Jo 

d6 cot[vr(iV+ \)e ie ] + f dB cot[-7r(iV + \) e w ] 



= 0. (B.9) 
Consequently, the complex integral on the RHS of (|B.7J) vanishes identically, q.e.d. 
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